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Abstract 

We produce an analog of the Davies-Unruh effect for superfluid helium four. There are two 
temperatures which result — one is associated with ordinary, or first, sound and the other with 
second sound. 
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Two of the most striking results of the investigations into quantum field theory on curved 

fa.- 



backgrounds are the Hawking effect [l| and the Davies-Unruh (D-U) effect 

la ISO 1 The 

Hawking effect predicts that black holes radiate particles thermally. The D-U effect predicts 
that a system interacting with some quantum field and following a particular non-inertial 
path undergoes transitions as if it were immersed in a thermal bath of excitations of the 
given quantum field, even though the quantum field itself is in a global Minkowski vacuum 
state. 

It is no coincidence that both results have thermal attributes. In each case that which 
is responsible for the thermal effects is an event horizon. 2 In the black hole example the 
horizon is a feature of the geometric structure of the underlying spacetime; for the Davies- 
Unruh scenario a horizon exists because of the particular trajectory of the system. The 
main effect of a horizon is to delimit regions of spacetime between which information cannot 
be exchanged. In 1973, Bekenstein p pointed out, for the particular example of black 
holes, that this inability to exchange information can bring about a violation of the second 
law of thermodynamics. He argued that this violation can be prevented if black holes have, 
in addition to their usual physical attributes, thermal characteristics (i.e. entropy). The 
very general nature of Bekenstein's reasoning suggests that it should always be the case that 
additional thermal properties result whenever an event horizon is present. 

In 1981, Unruh jj generalized the understanding behind these effects when he proposed 
an analog to the Hawking effect for irrotational perfect fluids. He considered an ordi- 
nary fluid which is steadily accelerated from zero velocity to supersonic flow. The place 
at which the transition to supersonic flow occurs, the sonic horizon, is the fluid analog of 
an event horizon. For example, the region of the fluid beyond the sonic horizon cannot 
communicate — through the exchange of fluid excitations (phonons) — with the starting place 
of the fluid flow (i.e. the fluid's "asymptotic infinity" where its velocity is zero). 3 Unruh 
found that the sonic horizon will emit phonons with a thermal distribution. The general- 



1 See also Fulling |j] for his preliminary contributions to the discovery of this effect. 

2 See Wald [5j for a suitable description of event horizons, in particular, the importance of causality (infor- 
mation cannot be exchanged instantaneously between two spatially separated points) in their definitions. 

3 Here, to "communicate" means to send information only through the exchange of phonons. Of course one 
could penetrate this horizon with something like a photon, say, but this should not ruin the event horizon 
nature of the sonic horizon for the problem at hand since it is interactions with the fluid excitations that 
dominate the end results. 
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ization occurs because here it is causality as defined by the paths that phonons take, not 
light, which underpins the physics. 

An important technical aspect of Unruh's derivation is that one can associate the fluid 
notions of sonic horizon and asymptotic infinity with an underlying "metric" in the math- 
ematics describing the system. Unruh manipulates this metric in precisely the same way 
that real metrics are manipulated (e.g., arbitrary coordinate transformations, construction 
of metric compatible covariant derivatives, etc.) to obtain sonic horizon thermal radiation. 
It is the main purpose of this work to apply Unruh's technique in order to demonstrate an 
analog of the D-U effect for superfluid He 4 . 

We use the two-fluid hydrodynamic model of superfluid He 4 developed by Landau 
The two fluids are the so-called normal- and super-fluids. The equations of motion for this 
model are (see Putterman j^|) 

= P,t + (pnVnj + PsVsj)j , (1) 
= (P S ),t + (P SV nj)j , (2) 
= V si)t + V sj V si)j + /i ; , (3) 

= (p n Vni + PsVsi) >t + (P$ij + PnVniVnj + PsVsiVsj)j , (4) 

where the repeated indices are summed over 1, 2, and 3, a subscript comma denotes partial 
differentiation, the subscripts n and s stand for normal- and super-fluid variables, respec- 
tively, any p, with a subscript or not, is a mass-density (p without a subscript is the total 
mass-density), Vj is the j th component of a three- velocity, p is pressure, s is the specific 
entropy, p is the chemical potential, and p s = p — p n . These variables are all functions 
of the Newtonian time t and the spatial coordinates Xj. The eight equations of motion in 
Eqs. (HJ)-(jlJ) are used to determine the eight independent degrees of freedom p, s, v n i, and 
v S i. To have a well-posed problem we must also supplement these equations with equations 
of state. Defining w; = v„j — and w 2 = WaW~, then the equations of state have the 
function depend! Q P = M = , M and , = M If T is the 
temperature, then the thermodynamic variables must also satisfy [2| 

dp = pdp + psdT + ^p n dw 2 . (5) 

The form of Eq. presupposes that the superfluid flow is irrotational, i.e. V x v s = 0. 
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This implies that there exists a scalar potential <p\ such that 



(6) 



As a consequence of Eq. (JBJ), and an assumption about initial conditions, there is another 
vector field Ai = (p n Wi)/(ps) which is irrotational [9|. Thus there is a scalar potential 02 
such that 

A = <h,i ■ (7) 

Because of the irrotational character of v S i and we really have only four independent 
degrees of freedom; therefore, some of the above equations are redundant. By inserting 
Eq. (jUJ) and Eq. (J7j) into Eqs. ([!])— (J3J) we can reduce the number of equations of motion from 
eight to four. Inserting Eq. © into Eq. Q yields 

= ($i,t + + l^j (8) 

so that 

= 0i,t + \Mi,3 +p + c 1 (t) , (9) 
where Ci(t) is an arbitrary function of t. Likewise, we can reduce Eq. (0J) to find 

= <j)2,t + —<faj<hj + (Pijfoj + T + c 2 (t) , (10) 

Pn 

where C2{t) is another arbitrary function of t. Now we have as our independent degrees of 
freedom (pi, <f) 2 , p, and s. Their equations of motion are Eqs. (JTJ) and (J2J) (with Eqs. (0) and 
© inserted appropriately), Eq. (j^J), and Eq. (fTUj). supplemented by the equations of state. 

To demonstrate our analog of the D-U effect we use a linearized form of the equations of 
motion about some background solution. The variables <f>i, 2 , P, s, p, T, p, and p n all take 
the form O = O + SO, where O represents a full solution, O represents the background, 
and SO is some variation away from the background. The prefactor S stands for a quantity 
which is "small" in the sense appropriate for this problem. We use a background which is 
homogeneous, isotropic, and static (i.e. all background functions <f>i, 02, etc. are independent 
of t and Xi). This implies that the thermodynamic variables p, T, p, and p n are functions 
of only p and s since the background velocities vanish. 

Before going any further it is convenient to relate some of the thermodynamic derivatives 
which appear during the linearization with more familiar quantities, i.e. the heat capacities 
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at constant volume Cy = Tds/dT\ p and constant pressure C p = Tds/dT\ p as well as the 
speeds of first and second sounds u\ = dp/dp\ s and u\ = (Ts 2 p s )/(Cyp n ), respectively l|. 
In terms of Cy, C p , ui, and u 2 we find 

dp 
ds 
dp 
dp 
dp 
ds 
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(13) 

(14) 



.Pn J \Pr, 

Since the background velocities vanish we do not require any specific formulas for the ther- 
modynamic derivatives of p n (these derivatives will not appear in the linearized equations). 

Inserting all variables which take the form O — O + SO into our equations of motion, 
keeping only terms linear in 5, and using the various thermodynamic quantities identified 
above, we find 
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We have done some further manipulations to arrive at Eqs. (|15 )) -([18 p : The linearized form 
of Eqs. Q and (fTUj) did not involve any derivatives of 5p or 5s so they could be solved in 
terms of derivatives of 5<pi and 5cf) 2 by purely algebraic means. 

We now show that some of the terms which appear in Eqs. (jl5j) - (fT£J) can be neglected 
for temperatures T = 10~ 3 K. In this temperature range the thermodynamics is dominated 
by phonons, the roton contributions being negligible Q, 11 ], and u\ = V3u 2 — 2.4 x 
10 4 cm/sec. Using the various equations of states for superfluid He 4 (see jlO, 11 1) and some 



basic thermodynamic identities, it can be shown that s = u\p n /T « 10 4 T 3 erg/(gm K 4 ) and 



(1 — Cy/Cp) 1 ' 2 ~ 10~ 2 T 2 K -2 . We have assumed that p s is constant (i.e. independent of 
T) so that dp/dT\ p « dp n /dT\ p . Using these approximations we find that the terms wf/uf, 
w\/{su^) ) and s[l — (w 2 , + w|)/m 2 ] are all much less than one for T < 10~ 3 K. Therefore, 
the equations for Scfti and 5(p2 decouple and take the following approximate form 

= — - S(/) a>tt + 8(j) adj , (19) 

where a = 1 or 2 is just a label, not a tensor index. 

Now we formally identify a metric from Eq. (|19|). This metric is obviously that associated 
with a flat spacetime: 

ds 2 = -u 2 a dt 2 + dx l dx i . (20) 

Note that time intervals are related to spatial intervals in this "spacetime" by the appropriate 
speed of sound u a . The important structure here is the sonic-cone given by the congruence 
of curves which satisfy ds 2 = 0. It is this structure which allows for the existence of a 
horizon in superfluid He 4 and underpins our Davies-Unruh analog. In contrast to real 
Minkowski spacetime, we see in this approximation to the dynamics of superfluid He 4 two 
sonic cones — one for first sound and another for second sound. 4 It is a general feature of 
Eqs. (fl3|) - (fT%|) that our system has two characteristic velocities for propagating excitations 
(it is not just a result of the approximations which generated Eq. (JUJ)) jo. Q liol 

To extract our analog of the D-U effect we will use the arguments developed by Unruh 
in Ref. j^. They involve the construction of "particle" detection devices. These devices 
register excitations of the quantum field under consideration, for us 5<p a . The detectors are 
to be accelerated, with constant proper acceleration a a , along the z-axis. (In this context, 
"proper" quantities are formed with the "metric" in Eq. (j2())l .) The worldlines traced out 
by the detectors have the following parametrization: 

z 2 - u 2 f = AR 2 , u a t = ztanh(r/2) , (21) 

and x = y = 0. The parameters r and R are the analogs of Rindler coordinates for the 
metric in Eq. (|20[) 



12j. In order to model how a detector registers excitations consider it to 



be a particle in a potential well. We will, for simplicity, let the potential be generated by 



4 We can anticipate similar results for third and fourth sound. Moreover, one can easily see that the 
Hawking effect analog will have more than one horizon (which will emit thermal excitations). 
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a box. The box must confine the particle and at the same time permit it to interact with 
the fluid. Let the interaction of the particle with 5(p a be via the potential 



where $ is the wave function of the particle and e a is some coupling constant. Note that 
our choice of Vi nt requires the particle to interact either with 8(pi or <50 2 but not both. 5 
Also, since the gradient of 5<f) a is a perturbation of a velocity, and velocities change sign 
under time reversal, Vi n t is time antisymmetric. This is not a serious problem here since our 
main result will be seen to be independent of the particular form of the interaction potential 
(i.e. no dependence on the coupling constant e a ). The overall behavior of the particle is 
governed by its Schrodinger wave equation 



R c of the box, a a = u 2 a /2R c is its proper acceleration, and the time derivative is with respect 
to the proper time of the particle (hence, the reason for the extra factor 2a a /u a ). The 
term proportional to m 2 just represents the force experienced by the particle because it is 
accelerating. 

Through some mechanism, which need not be discussed here, the detector is constrained 
to follow a R = R c = const trajectory in the < z < oo and — oo < t < oo region of a 
u a t — z spacetime diagram. As can be seen from Eq. (J2T|) this path is a hyperbola. The 
asymptotes of this hyperbolic path are generated by the sonic cone (i.e. the 45° lines in our 
spacetime diagram) going through the origin (t — z = 0). It is not very difficult to see that 
the "future" asymptote (i.e. the one with positive slope) does represent an event horizon 
since a fluid excitation produced "beyond" it can never be registered by our detector. 

The point of this construction is to see how the detector responds when it follows the 
given non-inertial path. A good indicator of this response is the transition probability per 
unit proper time P a j for the detector (i.e. "particle") to be excited from its ground state 

5 In supcrfluid He 4 it is possible to excite first sound modes without exciting second sound modes and vice 
versa. The reason for this is second sound modes are essentially thermal waves while first sound modes 
are associated with density disturbances. This follows from Eqs. (|15(1 and 1)16(1 since in the small T limit 
Sp oc 64>i.t and Ss oc 5(f>2,t- Also see Ref. Q. 
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Vint = e a 5<p a $ , 



(22) 




(23) 



where m is the mass of the particle, h ~ (R 2 — R 2 



)/R c is its proper distance from the center 



to some excited energy eigenstate with energy Ej. We will not present the details here, 
but rather we refer the reader to Unruh since the calculation of P a j follows exactly as 
presented there. 6 The final result is 

Pa] = exp(2iru a [Ej — E ] /a a )-l Fdet ' (24) 

where E is the ground state energy and F det just depends on the characteristics of the 
detector (e OJ etc.). Note that Eq. (j211) is precisely the form expected of a detector immersed 
in a thermal bath of excitations with temperature T a given by (restoring h and 

T a = -1 . (25) 
27T k b u a 

This is the desired result. (As claimed earlier, it is independent of e a .) For each speed 
of sound u a there is an associated temperature (cf. footnote 0]). Notice that this is the 
same formula obtained by Davies and Unruh except that the speed of sound u a replaces the 
speed of light. It is, up to factors of order one, what one would obtain from a dimensional 
analysis, using only the relevant physical parameters of this problem. Putting in the values 
for h, kb, and u a gives T a « 10~ 13 (a a / g) K where g ~ 10 3 cm/sec 2 . 

That we obtain two different temperatures is a necessary consequence of the fact that 
there are two different sonic-cone structures in superfluid He 4 . This is a required consistency 
check of the principle underlying our D-U analog — two sonic-cone structures necessarily 
implies two types of event horizons, which, in turn, suggests two different thermal effects. 

To arrive at Eq. (|23|) we have formally identified certain characteristics of superfluid 
He with an underlying metric present in the mathematics. We then formally used this 
metric in ways that are familiar from differential geometry and quantum field theory on 
curved backgrounds. At best, our derivation can only be considered heuristic since it is 
not obvious how to interpret the "proper time" ru a /2a a appearing in Eq. ([23)1 . It becomes 
especially problematic if one imagines a non-relativistic ion being the "detector" described 
by Eq. (J2EJ). The accepted understanding is that it is only Newtonian time, like t in Eqs. (JJJ- 
(J3J), which is the ion's proper time. If ru a /2a a has any real meaning it must be statistical, in 
essence, and perhaps, as Unruh states for his Hawking effect analog, "... (have) no physical 
significance except for the small-amplitude sonic motions of the medium" \jn. Whatever 



The important steps are those which produce Eqs. (3.1a) through (3.10) in Unruh 
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the interpretation, if one accepts the premise that horizons and thermal effects necessarily 
occur together, then one has little recourse but to believe that T a must have real, and not 
merely formal, significance. The problem, then, is to figure out how it can be detected. 
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Appendix: Note Added in 2005 

The work presented here is unpublished, and not otherwise available. The original 
preprint is of some historic interest: it was one of the earliest articles in the modern, post- 
1990 revival of "analog models" for curved spacetime. Every attempt has been made to 
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maintain the original form of the preprint. One typographical error, in Eq. (|1J). has been 
fixed. The referencing and footnoting are also different, since a modern typesetting program 
has been used. 
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